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ABSTRACT 

■ We present an analysis of optical spectroscopically-identified AGN down to a cluster 

magnitude of M* + 1 in a sample of 6 self-similar SDSS galaxy clusters at z 0.07. 
These clusters are specifically selected to lack significant substructure at bright limits 
in their central regions so that we are largely able to eliminate the local action of 
merging clusters on the frequency of AGN. We demonstrate that the AGN fraction 
increases significantly from the cluster centre to 1.5i?viriai, but tails off at larger radii. If 
only comparing the cluster core region to regions at ~ 2i?viiiai, no significant variation 
would be found. We compute the AGN fraction by mass and show that massive galaxies 
(log(stellar mass) > 10.7) are host to a systematically higher fraction of AGN than 
lower mass galaxies at all radii from the cluster centre. We attribute this deficit of 
AGN in the cluster centre to the changing mix of galaxy types with radius. We use 
the WHAN diagnostic to separate weak AGN from 'retired' galaxies in which the 
main ionization mechanism comes from old stellar populations. These retired AGN 
are found at all radii, while the mass effect is much more pronounced: we find that 
massive galaxies are more likely to be in the retired class. Further, we show that our 
AGN have no special position inside galaxy clusters - they are neither preferentially 
located in the infall regions, nor situated at local maxima of galaxy density as measured 
with E5. However, we find that the most powerful AGN (with [OHI] equivalent widths 
< —10 A) reside at significant velocity offsets in the cluster, and this brings our analysis 
into agreement with previous work on X-ray selected AGN. Our results suggest that if 
interactions with other galaxies are responsible for triggering AGN activity, the time- 
lag between trigger and AGN enhancement must be sufficiently long to obfuscate the 
encounter site and wipe out the local galaxy density signal. 

Key virords: galaxies: active — galaxies: evolution — galaxies: clusters: general 



1 INTRODUCTION 

Active galactic nuclei (AGN) are typically found inside mas- 
sive galaxies that exhibit significant, on-going or recent, star- 
formation (KaufFmann et al. 2003; Jahnke et al. 2004; Heck- 
man et al. 2005; von der Linden et al. 2010; Floyd et al. 
2012). The power source for AGN is expected to be gas ac- 
cretion on to a massive black hole (Lynden-Bell 1969) which 
suggests that black hole and galaxy spheriodal growth are 
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closely linked (cf. Richstone et al. 1998; KaufFmann et al. 
2003). 

Mergers have frequently been cited as a method to fuel 
AGN (e.g. Sanders et al. 1988) and a number of morpho- 
logical studies claim an excess of post-merger systems in 
their AGN samples (Bahcall et al. 1997; Canalize & Stock- 
ton 2001; Urrutia et al. 2008; Letawe et al. 2010; Smirnova 
et al. 2010). Given that the fuel source for AGN is in the 
gas phase, any physical mechanism that has the potential 
to disturb the morphology of a galaxy such as harassment 
(Moore et al. 1996) may also produce an enhancement of 
AGN activity - i.e. not simply mergers. Since such physical 
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Figure 1. A histogram of r-band magnitudes (with Poisson errors) for our sample is displayed in the left-hand panel. A line of best 
fit (dashed line) is fitted to the linearly increasing region of this plot (i.e. 16.0 < r < 17.0) which is used to create a completeness 
diagnostic plot (right-hand panel). The points in the right hand panel are the ratio of log(A'^) in the left-hand panel to line of best fit. 
The dotted lines enclosing these points denote the Ict error of the completeness values. At our adopted limiting magnitude of r = 17.77, 
the spectroscopy is still > 90 per cent complete. 



mechanisms can be tied to environment, AGN may there- 
fore be thought of as sign-posts to galaxy evolution in some 
circumstances (cf. Reichard et al. 2009). Indeed, there is ex- 
tensive literature supporting the idea that AGN (defined in 
various ways using different wavelengths) are influenced by 
environment. For instance, Kauffmann et al. (2004) report 
that the fraction of optical spectroscopic AGN is markedly 
different for galaxies in different density regimes (explicitly: 
the AGN fraction decreases as a function of increasing local 
galaxy density; see also Montero-Dorta et al. 2009). This is 
supported by a study of 51 galaxy clusters by Ruderman & 
Ebeling (2005) who find an excess of X-ray point sources 
within 3.5 Mpc of the centre of galaxy clusters in compar- 
ison to control samples which they attribute to AGN that 
have been triggered by close encounters with neighbouring 
galaxies. This broadly supports an increase of AGN fraction 
with increasing galaxy density, but they divide this excess 
into two regions: near the cluster core where galaxies are in- 
teracting with the brightest cluster member, and at around 
the virial radius whose excess they report is attributable to 
low energy collisions at the cluster-fleld boundary. Further, 
Popesso & Biviano (2006) detail an anti-correlation between 
cluster AGN fraction and cluster velocity dispersion (see also 
Sivakoff et al. 2008) . They contend that this anti-correlation 
indicates that the merger rate of clusters affects the AGN 
fraction since AGN are likely to have played a strong hand 
in heating the intra-cluster medium and thereby drive evo- 
lution in sub-groups that eventually form clusters (cf. Ma- 
mon 1992). From an investigation of the Abell 901/902 sys- 
tem, Gilmour et al. (2007) show that there is a deflcit of 
AGN in the highest density regions that supports the above 
works and other investigations of cluster AGN locations (see 
Gisler 1978; Dressier, Thompson & Shectman 1985; Cold- 
well, Martinez & Lambas 2002; Georgakakis et al. 2008; 
Gavazzi, Savorgnan & Fumagalli 2011; Pimbblet & Jensen 
2012). Naively, such a trend makes sense since galaxies in 
cluster centres would be more stripped of cold gas that can 



fuel an AGN than on cluster outskirts (see also Constantin 
et al. 2008; Lietzen et al. 2011). 

Equally, there is an increasing body of literature that 
indicates the opposite is true: environment plays little or 
no role in the frequency of AGN. Examining X-ray emis- 
sion from clusters. Miller et al. (2003) flnd no evidence for 
an enhanced cluster AGN fraction. This is supported by 
spectroscopic work on 8 clusters by Martini, Mulchaey & 
Kelson (2007) who detail the AGN fraction is no lower in 
cluster centres than a control fleld sample. Although some 
of this may be caused by mass selection effects (cf. Pasquali 
et al. 2009; Pimbblet & Jensen 2012), Haggard et al. (2010) 
demonstrate that there is no signiflcant difference in AGN 
fraction between cluster and field samples for a constrained 
range of absolute magnitudes. This is supported by von der 
Linden et al. (2010) who examine > 500 Sloan Digital Sky 
Survey (SDSS; Abazajian et al. 2009) clusters and find no 
trend in AGN fraction with distance from cluster centres 
(see also Atlee et al. 2011; Klesman & Sarajedini 2012). 

In this work, we present a new analysis of the AGN 
fraction's dependence on environment and mass to elucidate 
the issues summarized above using a sample of low redshift 
SDSS galaxy clusters that are free from known structure 
contamination. In Section 2 we detail the dataset that we 
use in this investigation. In Section 3, we compute how the 
AGN fraction varies with radius from the cluster centre and 
galaxy mass before discussing and summarizing our results 
in Section 4. Throughout this work, we adopt a standard, 
flat cosmology with Qm = 0.238, ~ 0.762 and Ho ~ 
73 km s"^ Mpc"^ (Spergel et al. 2007). 



2 DATASET 

Since mergers of galaxy clusters or sub-clusters may locally 
enhance AGN activity, we require a sample of clusters that 
are relatively free from such activity. In Pimbblet (2011), we 



© 0000 RAS, MNRAS 000, 000-000 



AGN fraction in galaxy dusters 3 



presented a sample of 14 such SDSS galaxy clusters derived 
from the earlier work of Plionis, Tovmassian & Andernach 
(2009) that have no merging or significant interaction with 
other (comparable) structures within the limits of SDSS ob- 
servations. In brief, Pliois et al. (2009) use the criteria and 
catalogue of Andernach et al. (2005) to generate a "clean" 
sample of Abell et al. (1989) clusters from. This consists of 
considering the velocity distribution of each cluster in turn 
and its flatness (Struble & Ftaclas 1994), and also removal 
of any cluster with > 1 X-ray peak. We refer the reader to 
Plionis et al. (2009) for a full description of this process. 

Here, we restrict the Pimbblet (2011) sample to 6 galaxy 
clusters that are within a narrow redshift slice (0.070 < z < 
0.084; Table 1). The reason for selecting such a narrow sam- 
ple to work with is to create a composite stacked cluster 
whose variation in absolute magnitude that corresponds to 
a given apparent magnitude completeness limit is small - no 
more than AMr = 0.4 (Pimbblet 2011). In terms of look- 
back time, the difference between our highest and lowest 
redshift clusters is ~ 0.1 Gyr. 

In-line with the SDSS spectroscopic limit (see Abaza- 
jian et al. 2009; Strauss et al. 2002), we use a limiting mag- 
nitude of r = 17.77 in this work. Fig. 1 demonstrates that at 
r = 17.77, SDSS spectroscopy is still > 90 per cent complete 
for our cluster sample (cf. Pimbblet & Jensen 2012; Jensen 
& Pimbblet 2012; Strauss et al. 2002). Further, Pimbblet 
(2011) notes that the sample covers little more than a factor 
of 2 in cluster mass. Combined with the small difference in 
look-back time, this ensures that the clusters in our sample 
are reasonably self-similar and are broadly at a comparable 
evolutionary stage. 

In Table 1, we give the global properties of the clusters 
used in this work, including mean recession velocity (cJ), 
cluster velocity dispersion (iTcz) and virial radius (-Rviriai). 
The former two are based on Miller et al. (2005) whilst the 
virial radius is computed from CTcz using the relation pre- 
sented by Girardi et al. (1998). Cluster membership is then 
simplistically defined to be all galaxies within ±3(Tcz of cz. 
The number of galaxies in each cluster within < Siiviriai is 
reported as N(< 37?viriai) in Table 1. These numbers are ap- 
proximately the same as those reported by Pimbblet et al. 
(2006; see their Table 2) for rich. X-ray luminous clusters at 
z - 0.1. 

To create our final sample, we stack all of our clusters 
together to form a composite sample. Analogous to Pimb- 
blet (2011), this is achieved by placing the clusters on to a 
common scale (i.e. i?viriai) and limiting our clusters to a com- 
mon absolute magnitude (i.e. the absolute magnitude corre- 
sponding to r- = 17.77 - Mr = -19.96 - at the redshift of 
our most distant cluster, Abell 1620; this corresponds to ap- 
proximately M* + 1 along the luminosity function according 
to the analysis of Jensen & Pimbblet 2012). As the terminal 
step, we select a mass limit for our clusters to prevent our 
sample being biased from having a low r-band limit that 
is coupled with a top-end mass limit (Holden et al. 2007; 
see also Pimbblet & Jensen 2012). To achieve this, we ex- 
amine plots of absolute magnitude versus stellar mass for 
each of our clusters and restrict our sample to those galax- 
ies brighter than A4r ~ —19.96 and more massive than the 
most massive galaxy log(stellar mass) — 10.4 at this limiting 
magnitude for our most distant cluster (Fig. 2). 

The final, bias-corrected composite sample consists of 
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Figure 2. Absolute r-band magnitude versus galaxy stellar mass 
(expressed as the logarithm of the solar mass of the galaxy) for 
the clusters in our sample. The vertical line denotes the spectro- 
scopic limit of r = 17.77 in our most distant cluster (Abell 1620) 
and the horizontal line denotes log(stellar mass) = 10.4 that cor- 
responds to the mass limit that we are ^complete to. These lines 
arc replicated for the other 5 clusters. We only use galaxies that 
are more massive and brighter than these limits for our subse- 
quent analysis to avoid biasing our sample. 



300 galaxies within i?viriai, and 701 galaxies within 3i?viriai 
from these 6 clusters. 



2.1 Sample Validation 

Before proceeding with our analysis, we elect to perform 
some validation tests on our sample to ensure that they are 
free of substructure - at least in their centres - as we have 
suggested. The need for such a validation step is our concern 
that the literature expresses both ambiguous and conflicting 
statements about some of the clusters in our sample. For 
example, Abell 1650 is noted by Einasto et al. (2012) as 
being a unimodal cluster, in agreement with Pimbblet et al. 
(2002) who term it a morphologically regular cluster using 
a large radii dataset. This is in contrast to Flin & Krywult 
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Table 1. The cluster sample used in this work. 



Name 


RA 


Dec 


cz 




^virial 


N(< 3KirM) 




(J2000) 


(J2000) 


(kms-l) 


(kms-1) 


(Mpc) 




A1205 


11 13 58.1 


+02 29 56 


22506 


938 


1.88 


110 


A1424 


11 57 26.4 


+05 05 52 


22764 


780 


1.56 


78 


A1620 


12 50 03.0 


-01 33 45 


25275 


1007 


2.01 


153 


A1650 


12 58 34.7 


-01 43 15 


25176 


864 


1.73 


117 


A1767 


13 36 31.6 


+59 08 51 


21111 


988 


1.98 


111 


A2670 


23 54 13.7 


-10 25 09 


22836 


976 


1.95 


132 



(2006) who note Abell 1650 has 'substructures in the cluster 

field'. 

Taking our bias-corrected sample, we apply the Dressier 
& Schectman (1988; DS) test for substructure to each clus- 
ter individually. The DS test is probably the most powerful 
available to detect substructure in arbitrary three dimen- 
sional datasets (Pinkney ct al. 1996) and we therefore con- 
sider it perfectly adequate for our sample validation check. 
Briefly, the DS test works by flnding the local mean velocity 
and standard deviation of the 10 nearest neighbours to a 
given galaxy and compares them to the global values for the 
cluster, such that: 

t-2 / A^local + 1 \ r/ \2 , / ^2^/1^ 

= I 2 I [(C^^local — C2:globalj +(,Criocal — f global) J(lj 

\ "■global / 

where the parameter S yields a measure of the deviancy 
of this sub-sample. A summed parameter of merit, A is then 
computed by summing all 5i terms in each cluster. To get a 
handle on the probability of A occurring, the velocity data 
are shuffled randomly between member galaxies 1000 times 
in a Monte Carlo fashion and the actual value of A is then 
compared to this ensemble. 

In Table 2, we give the values of P(A) for the clusters in 
our sample. All of our clusters are substructure free within 
IJ^viriai to the limits probed by our investigation. However, 
Abell 1620 stands out from the others as possessing signifi- 
cant substructure at high radii from the cluster centre. This 
agrees with the analysis of Burgett et al. (2004) who also 
suggested that the cluster may contain substructure from a 
complementary analysis of 2dF data. We analyze this clus- 
ter further in Appendix A. Despite the substructure at high 
cluster-centric radii, we retain Abell 1620 in our sample. 
We have experimented with removing this cluster from our 
sample and find that the effect on our primary results are 
negligible, but our uncertainties become fractionally larger. 

Finally, we issue the caveat that even though this sim- 
plistic test has indicated no substructure in the centres of our 
clusters, this does not preclude substructure at fainter mag- 
nitude limits arising from coherent, potentially low mass, 
infalling groups as might be expected (cf. Owers et al. 2011). 

2.2 AGN identification 

In order to identify which galaxies in our sample are AGN, 
we make use of a BPT diagram (Baldwin, Phillips &: Ter- 
levich 1981; see also Veilleux & Osterbrock 1987). The 
BPT plane consists of flux ratio of [NII]a6583/^?Q versus 
[OIU]\5007/Hfi. The measurements of equivalent widths for 
these lines are taken from Tremonti et al. (2004; see also 
www.mpa-garching.mpg.de/SDSS). In Fig. 3 we plot the po- 



Table 2. Results of the DS test using various cuts in radii. Results 
under 0.01 axe considered to indicate significant substructure. 



Cluster 


r < Rvirial 


P(A) 

r < 2Rvirial 


r < 3Rvirial 


A1205 


0.39 


0.12 


0.08 


A1424 


0.20 


0.03 


0.05 


A1620 


0.03 


< 0.01 


< 0.01 


A1650 


0.05 


0.05 


0.21 


A1767 


0.38 


0.34 


0.33 


A2670 


0.25 


0.31 


0.29 



sition of all galaxies in the composite cluster on the BPT 
plane that have S/N> 3 in the necessary lines. To differenti- 
ate AGN from galaxies that are simply star-forming, we use 
the demarcation curve of Kauffmann et al. (2003). The curve 
is a refinement of earlier work by Kewley et al. (2001) and 
yields 30 AGN within J?viriai and 81 AGN within 3i?viriai- 
We also define a composite sample - those galaxies that lie 
between the Kaufmann and Kewley curves - these galax- 
ies are weaker AGN whose host galaxies are star-forming. 
By implication, the use of the BPT diagram to select AGN 
means that our sample are composed of 'cold-mode' AGN 
(cf. Keres et al. 2005; Hopkins & Hernquist 2006). 

From the outset, we note that there is no significant 
optical colour difference (e.g. in g — r) between our AGN 
sample and the rest of the cluster population. This holds 
true even if we divide our sample by galaxy stellar mass. 



3 AGN FRACTION 

We now compute the cluster AGN fraction in two ways: 
by radius from the centre of the cluster and by galaxy 
mass. We note that the stellar masses of the galax- 
ies are from the SDSS value-added catalogue (www.mpa- 
garching.mpg.de/SDSS; see also Kauffmann et al. 2003). 

3.1 Fraction by Radius from Cluster Centre 

From the outset of our analysis, we note that we have se- 
lected the brightest cluster member as being the centre of 
our galaxy clusters. Although other choices could have been 
made, such as a luminosity-weighted centre or the peak X- 
ray flux location, we note that varying this choice does not 
alter the primary results presented in this work. We compute 
the AGN radial fraction in terms of RvMai (a more physi- 
cally meaningful scale than a fixed metric that uses Mpc; cf. 
Pimbblet et al. 2002) and plot the result in Fig. 4. The AGN 
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Figure 3. BPT plane for our composite cluster galaxy sample. All 
points have S/N> 3 in each line. The solid curve is the Kauffmann 
et al. (2003) demarcation line: galaxies above the curve are desig- 
nated AGN, those below are simply regular star-forming galaxies. 
Also shown is the Kewley et al. (2001) demarkation curve (dot- 
ted line). Galaxies between the two curves are composites: weaker 
AGN whose hosts are also star-forming. 



fraction is found to increase with distance from the cluster 
centre at a rate of d(fraction)/di?viriai = 0.018 ± 0.020. 

This rate of increase is not significant, which may ex- 
plain some of the differences reported in the literature con- 
cerning this fraction (for instance, contrast Gilmour et aL 
2007 and Kauffmann et al. 2004 with MiUer et al. 2003 and 
Martini et al. 2007). That said, if we restricted our analysis 
to the three innermost points of Fig. 4 (i.e. < 1.5-Rviriai) we 
find a significant gradient of d(fraction)/di?viriai ~ 0.064 ± 
0.021. For those studies that split the radial AGN fraction 
into bins similar to our analysis (Coldwell et al. 2002), sim- 
ilar results are obtained. 

Further, if we were to compare the AGN fraction within 
l.SiZviriai to that at 2.0 - 3.0 i^viriai, we find no significant 
difference (AGN fractions of 0.11 ± 0.02 and 0.11 ± 0.03, re- 
spectively). This may explain why investigations that that 
compare the 'cluster' environment to a control 'field' may be 
biased to finding no difference in AGN fraction. Moreover, 
the types of cluster used can also bias the measurement of 
AGN fraction - those investigations that use all types of 
cluster to probe AGN fraction could be biased by the pres- 
ence of substructure - a bias that the present work inten- 
tionally avoids - as could the use of very high central density 
clusters versus low density. 

The lines of best fit presented in Fig. 4 could also be 
over-simplifications of the situation. For example, Porter et 
al. (2008) report that there is an enhancement of specific 
star formation rate at ~few Mpc away from cluster centres 
caused by harassment inside galaxies being accreted along 
filaments of galaxies (see also Koyama et al. 2008; Perez et 



Figure 4. AGN fraction in the composite cluster as a func- 
tion of radius from the centre, with Poisson errors. Each bin 
covers 0.5-Rvirial except the last bin that covers li?virial to en- 
sure that each point in the plot has > 100 member galaxies (cf. 
Cameron 2011). A line of best fit to these points is overplotted as 
the solid line. This line has a gradient of d(fraction)/di?vii-iai = 
0.018 ± 0.020. We also display a line of best fit to the three 
points within 1.5-Rvirial (dotted line), which has a gradient of 
d(fraction)/(iRvirial = 0.064 ± 0.021. 

al. 2009). If true, then we may expect a similar enhancement 
of AGN fraction just beyond the virial radii of our clusters. 
Although Fig. 4 displays a local maxima in AGN fraction at 
1.75i?viriai, it is not significant - a larger sample of clusters 
and bona-fide filaments will be required to fully address this 
question. 

3.2 Fraction by Mass 

In Fig. 5 we compute the AGN fraction as a function of 
galaxy stellar mass for all galaxies within 37?viriai- We fit 
this data with a line of best fit and find that it has a signifi- 
cant gradient of d(fraction)/dlog(stellar mass) = 0.22 ±0.05 
- higher mass galaxies are significantly more likely to host 
AGN than their lower mass cousins (cf. Tanaka 2012; Pimb- 
blet & Jensen 2012; Xue et al. 2010; Brusa et al. 2009; Floyd 
et al. 2004; Dunlop et al. 2003; see also Best et al. 2005). 

To examine if a mass selection would affect the ra- 
dial cluster AGN fraction, we repeat our above radial anal- 
ysis for two mass bins in Fig. 6, split (arbitrarily) at 
log(stellar mass) = 10.7 to ensure that there are approx- 
imately equal numbers of galaxies above and below that 
mass in our sample. The more massive galaxies have a larger 
AGN fraction at all radii and there is a steady (but not 
very significant) increase in AGN from the cluster center to 
2-Rviriai bcfore it drops slightly lower again. Conversely, the 
lower mass galaxies do not vary in fraction significantly. The 
deficit of AGN in the centre of clusters may therefore sim- 
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Figure 5. As for Fig. 4, but as a function of galaxy stellar mass; 
each point has the same number of galaxies. The line of best 
fit has a gradient of c({fraction)/iilog(stellar mass) = 0.22 ± 0.05 
across the range studied. 



Figure 6. As for Fig. 4, but for galaxies with log(stellar mass) > 
10.7 (filled red circles, solid line) and < 10.7 (open blue triangles, 
dotted line). The higher mass regime contains a higher fraction 
of AGN at all radii. 



Table 3. Illustration of how the fraction of our most massive 
galaxies (log(stellar mass) > 11.0) in our sample changes as a 
function of radius. 



Radius 


Fraction 


(Rvirial) 




0-1 


0.19 ±0.02 


1-2 


0.16 ±0.03 


2-3 


0.14 ±0.03 



ply be a reflection of the changing mix of galaxy types (e.g. 
colour; morphology; mass) with cluster radius (cf. von der 
Linden et al. 2010). This is illustrated in Table 3, where we 
note the fraction of galaxies with log(stellar mass) > 11.0 in 
our sample steadily decreases with radius from the centre of 
our stacked cluster. This covariance of radius with mass is 
simply an expression of the well-known morphology-density 
relation (e.g., Dressier 1980; Dressier et al. 1997; Smith et 
al. 2005; see also Baldry et al. 2006). 

4 DISCUSSION 

4.1 AGN Phase-Space and Local Galaxy Density 

If AGN fractions are being enhanced in the cluster out- 
skirts by interactions with other galaxies, then an inves- 
tigation of their locations in {cz — c^)/(Jcz versus radius 
/ J?viriai phasB-space may reveal this. This phase-space is 
plotted in Fig. 7. We test if the AGN and other galaxy pop- 
ulations are distributed differently on this plane through a 
two-dimensional K-S test (Fasano & Franceschini 1987; Pea- 
cock 1983) and find that the two populations are the same 
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Figure 7. Phase-space diagram of AGN and other cluster mem- 
bers (dots) on a A(cz) versus radius plane. The AGN are coded 
according to their [OIII] emission. The largest blue points are 
AGN with an equivalent width of [OIII] < —10 A; the medium 
sized points have [OIII] equivalent widths in the range —1.5 - 
— 10 A; all the other (smallest) blue points have [OIII] equivalent 
width > -1.5 A. 
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Table 4. PrELction of retired AGN and galaxies from the total population taken from Fig. 8 as a function of radius and galaxy mass. The 
mass bins are chosen to have ap])roxiiiiatoly ocjiial j^alaxy iiiiiiibors in tho total y)oy)iilatioii. 



Bin 


Retired 
AGN fraction 


Retired fraction 
of total population 


0—1 Rvirial 


0.45 ±0.12 


0.37 ± 0.04 


1^2 ^virial 


0.34 ±0.10 


0.36 ± 0.04 


2—3 ^virial 


0.20 ±0.10 


0.41 ± 0.05 


log(stcllar mass) > 10.809 

10.575 < log(stcllar mass) < 10.809 

log(stellar mass) < 10.575 


0.51 ±0.10 
0.12 ±0.08 
0.12 ±0.08 


0.53 ± 0.05 
0.34 ± 0.04 
0.26 ± 0.03 



- i.e. the AGN are no more likely to be at large values of 
\{cz — cz)\/acz plus radius/_Rviriai (or conversely, centrally 
concentrated) than other cluster galaxies. This holds oven 
if we consider only galaxies more massive (or less massive) 
than log(stellar mass) = 10.7. Indeed, both the mean and 
median values of \{cz — cz)\/acz for the AGN and general 
cluster population are within la of each other; even within 
individual IRvMa.! radial bins. 

These results are in apparent disagreement with Haines 
et al. (2012; see especially their Fig. 6) who find that X- 
ray selected AGN are significantly more likely to be at the 
cluster infall regions than the general cluster population. 
Apart from using X-ray selected AGN, Haines et al. (2012) 
also use different cluster selection criteria: their clusters are 
more massive than the ones employed here and at higher 
rcdshift - the mean difference in lookback time between our 
sample and Haines et al. is « 1.7 Gyr - and the clusters may 
possess significant subclustering even at bright magnitudes. 
On the other hand, finding several AGN at comparatively 
low velocity offsets and radii is not absolute proof against 
them being an infalling population, as they can easily still 
appear at these locations (see Fig. 10 of Haines et al. 2012). 
An alternative hypothesis is that the AGN at low radii and 
velocity offsets arc 'retired galaxies' i.e. galaxies whose 
ionization mechanism is provided by old stellar populations 
(e.g. Cid Fernandes et al. 2010; 2011; Yan & Blanton 2012 
and references therein). In Fig. 8 we plot our sample using 
the WHAN diagnostic plot of Cid Fernandes et al. (2011). 
This plot is more 'efficient' than the standard BPT approach 
as it only uses two lines: Ha and [Nil]. Moreover, it is readily 
able to disentangle the so-called 'retired' galaxy population 
from weak AGN types. The WHAN diagram can also classify 
up to 50 per cent more of the emission line galaxy types than 
the BPT approach, and is therefore more able to distinguish 
LINERs from Seyferts, but no information from the BPT is 
'lost' in the move to the WHAN diagnostic. 

There are a number of noteworthy aspects of Fig. 8. 
At a basic level, it reflects already well-known results that 
actively star-forming galaxies reside at the outskirts of clus- 
ters (i.e. there are fractionally fewer circled dots in the star- 
forming corner compared the passive population), whereas 
the passive galaxies dominate the low velocity offset, low 
radii population (cf. Pimbblet et al. 2006; Pimbblet & Jensen 
2012 and references therein). Of our AGN population, some 
35 ± 6 per cent fall in to the 'retired' classification. These 
AGN are not preferentially situated at low radii and low 
velocity oflfsets: only 36 ± 10 per cent of the retired AGN 



satisfy such a criteria. This is, however, a slightly larger frac- 
tion than the other AGN classes: 14 ± 10 per cent of weak 
AGN and 16 ± 6 per cent of strong AGN reside at low radii 
plus low velocity offsets. These statistics are not significant 
enough to infer a duty cycle, but it is clear that all classes 
of AGN reside at all positions in our clusters. 

We interrogate Fig. 8 to determine how the retired frac- 
tion of both AGN and the entire galaxy population varies as 
a function of radius from the cluster centre and galaxy mass 
and present these results in Table 4. Whilst the retired frac- 
tion of all galaxies is approximately constant with radius, 
the retired fraction of AGN mildly decreases with distance 
away from the centre of our clusters. The trend with mass 
is steeper: high mass galaxies are much more likely to be 
in the retired class than lower masses. This result agrees 
with the contention of Lee et al. (2012) that the presence 
of a bar in a galaxy is no trigger for AGN activity (see also 
Combes 2003). Rather, the presence of bars and AGN is 
simply driven by host galaxy mass (e.g., Sheth et al. 2008; 
Nair & Abraham 2010). 

Taking Fig. 7 and Fig. 8 together, there is some sugges- 
tion that the disagreement with Haines et al. (2012) outlined 
above may not be so serious. By considering only the those 
AGN with equivalent widths of [OIII]< —10 A, we see that 
they seem to lie above \A{cz)\/acz ~ 1.2 (Fig. 7). This is 
reflected in Fig. 8: those AGN with the highest Ha emission 
values have larger velocity offsets (i.e. those points that are 
not circled in Fig. 8). This holds in all three panels of Fig. 8 
- i.e. at all radii. Haines et al. (2012) notes that their (pow- 
erful. X-ray) AGN reside on infalling caustic; hence if we 
only examine those powerful AGN in our sample, we come 
to an analogous conclusion. 

We can also examine if the AGN are preferentially in 
areas of high galaxy density by computing the local galaxy 
density for each galaxy in our sample. We choose E5 - the 
surface area on the sky that is occupied by a given galaxy to 
its tenth nearest neighbour - as our estimator of local galaxy 
density; this is effectively a probe of the internal densities 
of the dark matter halos (Muldrew et al. 2012). The plot of 
E5 as a function of radius from the cluster centre is shown 
in Fig. 9. There are a number of obvious AGN at high lo- 
cal galaxy density, but to test whether the AGN arc at a 
systematically higher value of S5 we compare bootstrapped 
mean Ss values of the AGN against the rest of the clus- 
ter members as a function of radius (Table 5) . This shows 
no significant difference between the two populations at any 
radii. 
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Figure 8. Diagnostic plot of Cid Fernandes et al. (2011) applied to our sample and split by radius. The larger blue points are our AGN, 
as defined by the Kauffmann demarkation line. The small dots are all our other galaxies, regardless of whether they have a high S/N ratio 
(i.e. unlike the definition of our AGN sample). Circled points denote low velocity offset galaxies with \ A{cz)\/acz < 0.5. Although ~35 
per cent of our AGN may be retired galaxies under this classification, those AGN with low radii and velocity offset are not preferentially 
retired. 



Therefore, if AGN are being triggered by encounters 
with other galaxies, they have since moved away from the 
site of the interaction suggesting a suitably long timescale 
between interaction and subsequent AGN enhancement. De- 
tailed computation of the value of this time-lag is beyond the 
scope of the present work, but we would expect it to be less 
than the time required for substructure to homogenize (e.g., 
Araya-Melo et al. 2009). Other studies (e.g. Schawinski et 
al. 2007; Shabala et al. 2012) have found that when inter- 
actions trigger AGN, it takes --100 - 200 Myr for an AGN 



to "switch on". Given a velocity dispersion of ~1000 km/s 
(cf. Table 1), each galaxy could move 100 - 200 kpc in this 
time. Therefore we regard it as comparatively easy to "wipe 
out" the local density enhancement signature by the time 
the AGN switches on. 

To round off this part of our discussion, we now cross- 
match our sample with the FIRST radio database to as- 
certain which of our galaxies would be classified as radio 
AGN (or 'hot-mode') and determine their radial fractions. 
We find 28 matches from FIRST using a cross-matching ra- 
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Table 5. Bootstrapped means and standard deviations of S5 
values for the AGN and other cluster members as a function of 
radius from the cluster centre. The two samples are statistically 
drawn from the same parent population. 



Radius AGN S5 Other S5 
(RvMal) (Mpc-2) (Mpc-2) 



0-0.5 


56.4 ±13.9 


39.1 ±5.6 


0.5-1.0 


18.4 ±2.7 


19.2 ± 1.5 


1.0-1.5 


35.2 ±20.8 


26.1 ±5.3 


1.5-2.0 


4.3 ±0.9 


5.4 ±0.6 


2.0-2.5 


8.8 ±3.2 


10.8 ±2.2 


2.5-3.0 


3.0 ± 1.5 


4.1 ± 1.7 



dius of 10 arcsec to our sample. To see whether the detected 
radio emission can arise due to just star formation, we exam- 
ine the SFR values for these galaxies from the value-added 
SDSS catalogues. We convert the reported SFR to a 1.4 
GHz radio luminosity using the formula of Yun et al. (2001) 
and determine how much higher the FIRST radio luminosity 
was than the value expected purely from star formation. Ra- 
dio AGN are then extracted as those galaxies whose FIRST 
radio luminosity exceeds that expected from the star for- 
mation within the SDSS fibre by a factor > la (where a is 
the standard deviation in the SFR estimate). This results 
in 15 radio AGN. This sample has a markedly high mass: 
the median is log(stellar mass/Af©) = 11.2 ± 0.3. The ra- 
dio AGN fraction for these galaxies goes as 0.013 ± 0.007, 
0.028±0.011, 0.027±0.012 for bins of li?viriai from the clus- 
ter centre. These small number statistics are hard to draw 
a meaningful conclusion from 



4.2 Composite AGN and AGN Power 

In Fig. 3, we identified not only the AGN (those galaxies 
above the Kaufmann demarkation line), but also a sample 
of composite weak AGN and star-forming galaxies (those 
between the Kaufmann and Kewley demarkation lines). In 
Fig. 10 we divide the radial AGN fraction in to the Kewley 
demarcated AGN, and the composite sample to ascertain 
if the composite sample are driving any of the trends seen 
above. The composite sample appears to follow quite a flat 
distribution. Meanwhile, the Kewley et al. (2001) defined 
AGN show a steeper initial variation with increasing radius 
which plateaus quickly. Both samples are consistent with 
the trend depicted in Fig. 4 for the AGN fraction gradient. 
However, if we consider only the inner points (i.e. the dotted 
line in Fig. 4), we see that this is more consistent with the 
Kewley defined AGN. 

If there is a bona-fide radial trend of AGN fraction with 
radius, we may be able to see this reflected in the [OIII] line 
strength which is a proxy for AGN power. In Table 6, we list 
the median equivalent widths of [OIII] for our original AGN 
(i.e. Kauffmann delineated) with radius (see also the differ- 
ent point sizes in Fig. 7 and Fig. 9). For comparison, we also 
detail the same measurement for all galaxies in our sample. 
No significant radial trend is observed for either sample. In 
the case of the AGN, this may simply be because we lack the 
numbers to detect such a trend since we only have 84 AGN 
within 3-Rviriai (and a maximum of 32 galaxies in a IMpc 



radius / R„, 



Figure 9. Local galaxy density, S5 as a function of radius from 
the cluster centre. The large blue points denote the AGN, coded 
as per Fig. 7; small dots are the other galaxies. The AGN sta- 
tistically occupy the same region of this parameter space as the 
other cluster members (Table 5). 




Radius / R^i, 



Figure 10. As for Fig. 4, but split on the basis of whether the 
galaxies falls in the 'composite' AGN+star-forming sample (black 
open triangles), or the Kewley et al. (2001) definition of AGN 
(green filled circles). The two samples are fractionally offset in 
radius from each other for clarity. 
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Table 6. Variation of [OIII] equivalent width (EW) with radius. 
The quoted uncertainties are bootstrapped standard deviations. 



Radius 


AGN Sample 


All Galaxies 


(^virial) 


Median [OIII] EW 


Median [OIII] EW 




(A) 


(A) 


0-1 


-1.3 ±0.2 


-0.4 ±0.1 


1-2 


-1.6 ±0.7 


-0.5 ±0.1 


2-3 


-1.6 ±0.7 


-0.6 ±0.1 



bin in Table 6). A larger sample will be required to inves- 
tigate this. Wc can, however, infer that we arc not missing 
any AGN because they are simply too weak and therefore 
below the BPT detection threshold. Hence the change in 
AGN fraction is a bona-fide change in the duty cycle rather 
than AGN luminosity. 



4.3 AGN Colour 

If AGNs axe preferentially associated with bright galaxies 

that are presumably morphologically late-type, then deter- 
mining the fraction of AGNs in radial bins in our sample 
should mimic the radial distribution of late-type galaxies 
in clusters - i.e. the morphology-density relation (Dressier 
1980; see above) and therefore not address whether a spe- 
cific mechanism is tied to AGN activation or quenching. To 
partially resolve this, we could attempt to morphologically 
classify all of our sample through using (e.g.) GalaxyZoo 
(Lintott et al. 2008) in an analogous way to Pimbblet & 
Jensen (2012). But this approaches contain problems: prin- 
cipal among them being a large fraction of 'uncertain' clas- 
sifications that could skew an analysis. To combat this, we 
construct a colour-mass diagram for our sample, divided by 
radius to the cluster centre (Fig. 11). We divide this dia- 
gram up in to two halves: a red sequence and a blue cloud 
component. This is done by eye, choosing a lino that divides 
the two reasonably cleanly. Despite the arbitrary nature of 
this approach, it serves our purpose of creating two cate- 
gories of galaxies (i.e. early and late types) and we note that 
the gradient is consistent with the colour-magnitude rela- 
tions presented in earlier works (Pimbblet et al. 2006; 2002) 
and the lower envelope limit of such fitted colour-magnitude 
relations. However, we caution that we have made no at- 
tempt to correct the (g-r) colours in Fig. 11 for AGN blue- 
ing: given the use of SDSS model magnitudes, we suggest 
that this ofi^cct would be small and only likely to affect the 
strong (i.e. bluest) AGN population since the weak and re- 
tired AGN are already predominantly residing on the red 
sequence (Fig. 11). 

From Fig. 11, we see that any enhancement of the red 
sequence AGN fraction at low radii appears to be purely 
driven by massive galaxies. We quantify this in Table 7 
where we detail the AGN fractions not only above and below 
the red sequence envelope cut-off, but also divided by mass. 
Although the uncertainties on these numbers are large (too 
large to infer statistically significant trends) , by considering 
galaxies lying within and outside 1 i?viriai of the cluster cen- 
tre, it is tempting to speculate that red massive AGN and 
blue low-mass AGN may have a radial dependence whereas 



the red sequence low mass AGN may not. If so, this may im- 
ply a common AGN triggering mechanism such as gas-rich 
interaction. Hence if a low-mass galaxy underwent such an 
interaction, it would necessarily become blue due to the par- 
allel star formation. This would not be the case for a massive 
galaxy since it will have a lower specific star formation rate 
which in turn would correlate with the (g-r) colour. A larger 
sample of clusters is required to unambiguously address this 
issue. 

Finally, we note that disregarding the division by mass, 

at all radii the fraction of AGN in the blue category is twice 
that in the red. This is broadly consistent with our earlier 
results that galaxies which show signatures of recent inter- 
actions show elevated levels of both AGN activity and blue 
colours (Shabala et al. 2012; Kaviraj et al. 2012). 



5 SUMMARY 

In summary, we have investigated the AGN fraction in 6 

'clean' galaxy clusters down to « M* + 1 as a function of 
both mass and radius from the cluster centre. Our main re- 
sults are: 

(i) The radial AGN fraction increases steeply in the central 
1.5 /?viriai of the compositc cluster, but flattens off quickly 
and even decreases beyond this radius. If one were to com- 
pare the central regions of clusters with field samples, then 
no difference would be found on the basis of this work. 

(ii) The AGN fraction by mass shows a significant trend such 
that more massive galaxies arc more likely to host AGN. In- 
deed, massive galaxies host more AGN at all radii from the 
cluster centre. The reported deficit of AGN in cluster cen- 
tres may therefore simply be a product of the changing mix 
of galaxy types with radius. 

(iii) Retired AGN (found using the WHAN diagnostic) are 
found at all radii in the cluster, but their response to mass 
is much more pronounced: we find that massive galaxies are 
more likely to be in the retired class. 

(iv) AGN have no preferential position inside galaxy clus- 
ters (either with regard to infalling status, or enhanced local 
galaxy density). This conclusion can be brought in to line 
with studies of X-ray AGN (e.g. Haines et al. 2012) by con- 
sidering only the most powerful optical AGN. These galax- 
ies avoid \ A{cz)\/acz < 1.2 and may therefore reside on the 
cluster caustics (infall regions) as demonstrated in Haines et 
al. (2012). 

(v) If interactions with other galaxies trigger AGN activity, 
then the time-lag between the trigger and AGN enhance- 
ment must be sufficiently long to mask the site of the en- 
counter and eliminate any signal in local galaxy density. 

Our favoured scenario for AGN triggering remains a 
gas-rich interaction, although increased numbers of galaxies 
are required to produce better statistics to firm this specu- 
lation up. 
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to their WHAN classification from Fig. 8 (crosses for strong AGN; filled triangles for weak AGN; filled circles for the retired class). The 
dashed diagonal line is our approximation for the division between red sequence (above the line) and blue cloud galaxies (below the line). 
The numbers show the corresponding AGN fractions above and below this line. These fractions remain statistically constant with radius. 



S.S.S. thanks the Australian Research Council for the 
award of a Super Science Fellowship at the University of Tas- 
mania. D.J.E.F. acknowledges support from the Australian 
Research Council via Discovery Project grant DP110102174. 

Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating In- 
stitutions, the National Science Foundation, the U.S. De- 
partment of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 



Planck Society, and the Higher Education Funding Council 
for England. The SDSS Web Site is http://www.sdss.org/. 

The SDSS is managed by the Astrophysical Research 
Consortium for the Participating Institutions. The Partic- 
ipating Institutions are the American Museum of Natu- 
ral History, Astrophysical Institute Potsdam, University of 
Basel, University of Cambridge, Case Western Reserve Uni- 
versity, University of Chicago, Drexel University, Fermilab, 
the Institute for Advanced Study, the Japan Participation 
Group, Johns Hopkins University, the Joint Institute for 



© 0000 RAS, MNRAS 000, 000-000 



12 K. A. Pimbblet 



Table 7. AGN fractions derived from Fig. 11, divided by radius, colour and mass. The colour dividing line for blue cloud and red 
sequence is taken as the dashed line from Fig. 11, whereas the division in mass between massive and low mass galaxies is at log(stellaj: 
mass)=11.0. 



Radius 


Massive 


Low Mass 


Massive 


Low Mass 


(/?virial) 


Red Sequence 


Red Sequence 


Blue Cloud 


Blue Cloud 


0-1 


0.20 ± 0.06 (10/51) 


0.06 ±0.02 (13/207) 


0.00 ±0.00 (0/2) 


0.15 ±0.06 (6/40) 


1-2 


0.28 ± 0.09 (10/35) 


0.08 ±0.02 (10/128) 


1.00 ±1.00 (1/1) 


0.22 ±0.07 (11/50) 


2-3 


0.19 ±0.11 (3/16) 


0.05 ± 0.02 (6/120) 


0.00 ± 0.00 (0/4) 


0.23 ±0.07 (11/47) 
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Figure 12. Space density of our bias-corrected galaxy sample for Abell 1620. The galaxies (white crosses) have been smoothed with an 
(arbitrary) Gaussian kernel of 0.1 deg in length and the density represented by the variation in colour (see the side colour bar). The centre 
of the plot corresponds to the position of the brightest cluster galaxy in Abell 1620, and is also the luminosity-weighted centre. Two 
other peaks in the spatial distribution are seen at 12h48m02.7s -01d39ml0s and 12h48m57.7s -01d55m33s that we identify as SDSS-C4 
1010 and NSC J124857-015532 respectively. Although these secondary peaks appear significant in space density, their masses are less 
than the primary A1620 cluster, in agreement with Burgett et al. (2004). 
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APPENDIX A: ABELL 1620 

Abell 1620 is noted in Section 2 as possessing significant 
substructure at > 2J?viriai, in agreement with Burgett et al. 
(2004). 

In Fig. 12 we plot a smoothed surface density of galaxy 
members for this cluster. This figure reveals two overdensi- 
ties of galaxies to the South East of Abell 1620 proper. We 
identify these two overdensities as SDSS-C4 1010 (Miller et 
al. 2005) at 12h48m02.7s -01d39ml0s, and NSC J124857- 
015532 (Gal et al. 2003) at 12h48m57.7s -01d55m33s using 
NED. 

These sub-clusters (i.e. groups) are notable in our anal- 
ysis. In Fig. 9, there is a local peak in S5 at 1.2_Rviriai- 
1.5i?viriai. Wc associate this peak with these two groups. Of 
note, there are two AGN contained in this peak (i.e. with 
log(E5) > 1.8). It is these AGN that result in the enhanced 
E5 average value noted in Table 5 at these radii. Wc ex- 
plicitly note that the mass of these secondary peaks arc less 
than the primary A1620 peak (as confirmed by Burgett et al. 
2004). We further emphasize that removal of Abell 1620 from 
our analysis does not change the primary results contained 
in our work. Therefore, wc regard these two overdensities to 
have had negligible effect on the AGN contained within it 
(i.e. no enhancement), in-line with our conclusions. 
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